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ABSTRACT 


We aim to gain insight into the effect of network and faculae on solar irradiance from their apparent intensity. Taking full-disc 
observations from the Solar Dynamics Observatory, we examined the intensity contrast of network and faculae in the continuum 
and core of the Fe I 6173 À line and 1700 À, including the variation with magnetic flux density, distance from disc centre, nearby 
magnetic fields, and time. The brightness of network and faculae is believed to be suppressed by nearby magnetic fields from its effect 
on convection. The difference in intensity contrast between the quiet-Sun network and active region faculae, noted by various studies, 
arises because active regions are more magnetically crowded and is not due to any fundamental physical differences between network 
and faculae. These results highlight that solar irradiance models need to include the effect of nearby magnetic fields on network and 
faculae brightness. We found evidence that suggests that departures from local thermal equilibrium (LTE) might have limited effect 
on intensity contrast. This could explain why solar irradiance models that are based on the intensity contrast of solar surface magnetic 
features calculated assuming LTE reproduce the observed spectral variability even where the LTE assumption breaks down. Certain 
models of solar irradiance employ chromospheric indices as direct indications of the effect of network and faculae on solar irradiance. 
Based on past studies of the Ca II K line and on the intensity contrast measurements derived here, we show that the fluctuations in 
chromospheric emission from network and faculae are a reasonable estimate of the emission fluctuations in the middle photosphere, 
but not of those in the lower photosphere. The data set, which extends from 2010 to 2018, indicates that intensity contrast was stable 


to about 3% in this period. 
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1. Introduction 


The variation in the solar radiation that enters the Earth's at- 
mosphere is a key consideration in the climate change debate 
(Gray et al. 2010). For this reason, climate simulations require 
historical solar irradiance, defined as the Earthward solar radia- 
tive flux at one AU, as input. This is provided by models of so- 
lar irradiance variability because measurements only date back 
to 1978, and even then have gaps and gross instrumental issues 
(Ermolli et al. 2013; Kopp 2014). The variation in solar irradi- 
ance at solar rotation to cycle timescales is believed to be domi- 
nantly driven by photospheric magnetism (Domingo et al. 2009; 
Yeo et al. 2017b). Specifically, it is thought to be the outcome 
of the changing prevalence and distribution of dark sunspots and 
bright network and faculae on the solar surface. Models based 
on this assumption have been successful in replicating most of 
the features in measured solar irradiance (Solanki et al. 2013). 
For example, the SATIRE-S model (Wenzler et al. 2006; Ball 
et al. 2012; Yeo et al. 2014) reproduces over 90% of the observed 
variability in the wavelength-integrated or total solar irradiance 
since 1978. Of course, a firm understanding of the radiant prop- 
erties of sunspots, network, and faculae is pertinent to the proper 
modelling of solar irradiance variability. While sunspot radiance 
is well constrained (Maltby et al. 1986; Collados et al. 1994; 
Unruh et al. 1999), the same cannot be said of network and fac- 
ulae, mainly because of the challenge in observing such small- 
scale structures. 


The radiant properties of network and faculae are often 
probed by taking simultaneous solar magnetograms and inten- 


sity images, and studying the variation in network and facular in- 
tensity contrast with distance from disc centre and magnetogram 
signal (Topka et al. 1992, 1997; Lawrence et al. 1993; Ortiz et al. 
2002; Yeo et al. 2013). As magnetic flux tubes tend towards 
a surface normal orientation as a result of magnetic buoyancy, 
the viewing geometry is effectively a function of distance from 
disc centre. Above a certain minimum, magnetic flux tubes ex- 
hibit similar kilogauss magnetic flux densities regardless of size 
(Solanki & Schmidt 1993; Solanki et al. 1999). This means that 
the magnetogram signal at network and faculae is indicative of 
the magnetic filling factor. Moreover, the magnetic flux tube size 
rises on average with magnetogram signal (Ortiz et al. 2002; Yeo 
et al. 2013). To study the intensity contrast of network and facu- 
lae as a function of distance from disc centre and magnetogram 
signal is to study how their radiance varies with viewing geom- 
etry, magnetic filling factor, and flux tube size. 


The Topka et al. (1992, 1997) and Lawrence et al. (1993) in- 
vestigations were based on scans of various disc positions made 
with the Swedish Vacuum Solar Telescope (SST, Scharmer et al. 
2003). The scans varied in spatial resolution (0.3 to 0.6 arcsec) 
and wavelength, making it difficult to constrain the disc centre- 
to-limb variation (CLV) in intensity contrast. In addition, taken 
from the ground, SST observations suffer from atmospheric see- 
ing. Ortiz et al. (2002) used full-disc observations from the 
Michelson Doppler Imager on board the Solar and Heliospheric 
Observatory (SoHO/MDI, Scherrer et al. 1995). These obser- 
vations are not only seeing-free because they are taken from 
space, but cover the entire solar disc at a uniform if relatively 
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poor spatial resolution (4 arcsec). This allowed the authors to 
derive an empirical relationship describing the variation in net- 
work and facular intensity contrast with distance from disc cen- 
tre and magnetogram signal. In the various studies, the authors 
examined the intensity contrast in the nearby continuum of pho- 
tospheric lines. The spectral line intensity contrast is also of in- 
terest, not least because studies have found evidence that both 
continuum and spectral line changes might be relevant to solar 
irradiance variability (Mitchell & Livingston 1991; Unruh et al. 
1999; Preminger et al. 2002). However, early investigations into 
the line core intensity contrast of network and faculae were lim- 
ited by insufficient data, such that they were unable to constrain 
either the CLV or the magnetogram signal dependence with 
confidence (Frazier 1971; Walton 1987; Lawrence et al. 1991; 
Title et al. 1992). Taking observations from the Helioseismic 
and Magnetic Imager on board the Solar Dynamics Observatory 
(SDO/HMI, Scherrer et al. 2012), Yeo et al. (2013), abbreviated 
here as Yeal3, examined the intensity contrast of network and 
faculae in both the continuum and the core of the Fe I 6173 A 
line. HMI returns seeing-free full-disc observations of the mag- 
netic flux density and the continuum and line core intensity at 
a uniform, intermediate spatial resolution (0.9 arcsec), allowing 
Yea13 to circumvent the limitations of the earlier studies. Using 
a simple model based on the observed CLV and magnetograms 
signal dependence of the continuum and line core intensity con- 
trast, the authors confirmed that solar irradiance variability is in- 
deed the sum manifest of continuum and spectral line changes. 


Comparing network and facular intensity contrast by looking 
at quiet and active regions separately, Lawrence et al. (1993) and 
Kobel et al. (2011) noted that above a certain minimum magne- 
togram signal level, the quiet-Sun network is brighter than active 
region faculae at a similar magnetogram signal level. Ortiz et al. 
(2002) and Yea13 found that the intensity contrast per unit mag- 
netic flux density decreases with increasing magnetic flux den- 
sity. Taking into consideration that in their analyses, the points 
with lower magnetic flux densities correspond mainly to network 
and higher levels to faculae, both studies concluded that network 
is, on a per unit magnetic flux density basis, brighter than facu- 
lae. The various investigations point to network having a greater 
heating efficiency than faculae. This phenomenon is well sup- 
ported by observations (Solanki & Stenflo 1984; Solanki 1986; 
Keller et al. 1990) and magnetohydrodynamics (MHD) simula- 
tions (Vógler 2005; Criscuoli 2013). This has been attributed 
to the assumption that facular regions are more magnetically 
crowded; the more magnetically crowded a given area, the more 
convection is impeded and the less efficient the radiative heat- 
ing of magnetic flux tubes by surrounding convection (Ishikawa 
et al. 2007; Kobel et al. 2012; Romano et al. 2012; Criscuoli 
2013). Significantly, by studying areas with varying degrees of 
magnetic activity, Kobel et al. (2012) found that the intensity 
contrast at a given magnetogram signal level and the Doppler ve- 
locity dispersion, which is an indication of convective strength, 
both decline gradually with the area-averaged magnetic flux den- 
sity. 

Here, we present a follow-up study to Yea13, who took HMI 
observations and examined the CLV and magnetogram signal 
dependence of network and facular intensity contrast in the con- 
tinuum and core of the Fe I 6173 A line. Making use of similar 
HMI data and simultaneous observations from the Atmospheric 
Imaging Assembly (AIA, Lemen et al. 2012), also on board 
SDO, we study the intensity contrast in the AIA 1700 À channel 
along that in the Fe I 6173 À line. In another advance, we exam- 
ine not only the CLV and magnetogram signal dependence, but 


also the variation in intensity contrast with the average magnetic 
flux density over the neighbourhood of each network and facu- 
lar image pixel, called here the local magnetic flux density, and 
time. The objectives of the study are listed below. 


— We aim to advance what Yeal3 understood about network 
and facular radiance from the intensity contrast in the Fe I 
6173 À line by a comparison to the contrast in the AIA 1700 
À channel. 

— Kobel et al. (2012) studied the disc centre. We extend this 
earlier study by examining the variation in intensity contrast 
with local magnetic flux density across the entire solar disc. 
This is to shed light on how the observation that network has 
a higher heating efficiency than faculae might be appropri- 
ately described in solar irradiance models. 

— HMI and AIA have been in operation since 2010, almost the 
length of an entire solar cycle. We investigate if network and 
facular radiance might vary with cycle phase by looking at 
intensity contrast as a function of time. Clearly, this also has 
implications for the modelling of the effect of network and 
faculae on solar irradiance. 

— Certain solar irradiance models adopt chromospheric indices 
as direct indications of the effect of network and faculae on 
solar irradiance (e.g. Morrill et al. 2011; Thuillier et al. 2012; 
Coddington et al. 2016; Yeo et al. 20172). Taking the inten- 
sity contrast measurements derived here and what is known 
about Ca II K radiance into consideration, we discuss the 
validity of this assumption. 


In the following section, we describe the HMI and AIA observa- 
tions we employed. Then, we detail how we isolate the network 
and faculae in the data set (Sect. 3.1), and derive the correspond- 
ing intensity contrast (Sect. 3.2) and local magnetic flux density 
(Sect. 3.3). We present and discuss the results in Sect. 4, before 
giving a summary of the study in Sect. 5. 


2. Data 
2.1. HMI 


The HMI (Scherrer et al. 2012) records full-disc polarised fil- 
tergrams at the Fe I 6173 À line at 1.875s intervals on two 
4096 x 4096 pixel CCDs. The pixel scale and spatial resolution 
is 0.5 and 0.9 arcsec, respectively. A range of data products is 
generated from the filtergram sequence, including 45s and 720s 
cadence Dopplergrams, longitudinal magnetograms, and maps 
of the continuum intensity, and the depth and width of the Fe I 
6173 À line. The line core intensity is given by the difference 
between the continuum intensity and line depth maps. Vector 
magnetograms are also generated at 720s cadence. 

For this study, we made use of 45s longitudinal magne- 
tograms, continuum, and line core intensity images. The contin- 
uum and core of the Fe I 6173 A line is formed in the lower pho- 
tosphere (around 20 km, Norton et al. 2006) and over the middle 
photosphere (the contribution function is centred on 250 km with 
a full width at half-maximum, FWHM, of 300 km, Fleck et al. 
2011), respectively. 

The longitudinal magnetogram gives the mean line-of-sight 
magnetic flux density over each resolution element, Bios. As 
magnetic flux tubes tend towards a surface normal orientation, 
we can approximate the pixel-averaged magnetic flux density, 
B by the quotient Bios/u, where u is the cosine of the heliocen- 
tric angle. The 1c noise level of HMI 45s longitudinal magne- 
tograms and 720s vector magnetograms is about 10 G (Liu et al. 
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2012) and 100 G (Hoeksema et al. 2014). The elevated noise 
level of HMI vector magnetograms, arising from the uncertainty 
in measured Stokes Q and U, is the reason why we choose to 
estimate B from Bios over taking it directly from vector magne- 
tograms. 


2.2. AIA 


The AIA (Lemen et al. 2012) comprises four 4096 x 4096 pixel 
CCD telescopes that image the Sun in ten wavelength channels. 
The pixel scale and spatial resolution are 0.6 and 1.5 arcsec. 
Seven channels are centred on extreme ultraviolet (EUV) lines 
that are formed in the transition region and corona. The other 
three channels are at nominally 1600 Á, 1700 A, and 4500 A. 
Full-disc images are recorded in each EUV channel at 12s inter- 
vals, the 1600 A and 1700 A channels at 24s intervals, and the 
4500 A channel hourly. 

Here, we examine the intensity contrast of network and fac- 
ulae in the 1700 A channel. The 1700 A channel captures con- 
tinuum emission, mainly from neutral silicon bound-free transi- 
tions, that forms over the middle photosphere (the response func- 
tion is centred on 360 km with an FWHM of 325 km, Fossum & 
Carlsson 2005). 

We excluded the other AIA channels from this study for the 
following reasons. The EUV channels are formed in the transi- 
tion region and corona, where magnetic structures are no longer 
co-spatial to their photospheric footpoints, and are therefore not 
suitable for such a study. The 1600 A channel captures contin- 
uum emission that forms in the upper photosphere and C IV 
emission that forms in the transition region (Handy et al. 1999; 
Lemen et al. 2012), which complicates the interpretation of the 
apparent intensities. AIA 4500 A images have severe artefacts 
that are caused by damage to a filter.! 


2.3. Data selection and reduction 


HMI and AIA have been in continuous operation since April 
2010. We made use of data from 100 days between June 2010 
and October 2018, selecting from each calendar month the day 
with the highest sunspot number (Clette et al. 2016). 

For each data day, we identified a 315s period over which 
HMI and AIA operated without interruption. We took the seven 
HMI 45s longitudinal magnetograms and continuum and line 
core intensity images from this interval, and the 11 AIA 1700 
A images recorded over the same period. We rotated each mag- 
netogram to the middle of the 315s period and took the average, 
likewise for the other data products. The result is a 315s lon- 
gitudinal magnetogram and intensity image at each passband, 
which are co-temporal (Fig. 1). This averaging suppresses the 
variability from noise, p-mode oscillations (Krijger et al. 2001), 
and acoustic waves excited in the convection zone (Wedemeyer 
et al. 2004). 

We co-registered the HMI and AIA observations using the 
aia_prep.pro routine (SolarSoft IDL). This resamples the HMI 
data to the AIA image pixel scale. It is known that the co- 
registration by this routine is only accurate to about | arcsec be- 
cause of short-term variability in the pointing of AIA and HMI, 
possibly from thermal flexing, which is not captured in the data 
header (see Brice Orange et al. 2014, and references therein). We 
refined the alignment between the HMI and AIA data by com- 


! Reported by Dean Parnell on the official SDO mission weblog, 
sdoisgo.blogspot.de, in the entry dated 17 December 2013. 


paring the absolute magnetogram signal to the 1700 A intensity 
at disc centre, taking care to exclude sunspots, and finding the 
shift to the magnetogram that maximises the cross-correlation 
between the two. The required shift varied between 0.3 and 1.2 
arcsec over the data days. 


3. Analysis 
3.1. Identification of network and faculae 


Network and faculae were identified following Yea13. First, we 
determined the noise level of the HMI 315s longitudinal mag- 
netograms, ca, as a function of position on the solar disc by 
the exact procedure described in Yeal3. (We cannot employ the 
noise level estimate from this earlier study here as it was de- 
termined for HMI observations at the original spatial resolu- 
tion.) The noise level increases radially, from an average of 4.5 
G around disc centre (0.9 < u < 1.0) to 6.4 G near the limb 
(0.1 < u < 02). 

For each data day, we took the points on the solar disc where 
Bios > 30B, in the longitudinal magnetogram to harbour mag- 
netic activity. The weakest magnetic flux density detectable, Bin; 
is given by 30g, /u. This rises from 14 G around disc centre to 
127 G near the limb from the combination of the disc centre- 
to-limb increase in the noise level and foreshortening. The con- 
tinuum intensity image is normalised by the quiet-Sun level and 
flat-fielded (detailed next in Sect. 3.2). The active pixels where 
the normalised continuum intensity is below 0.89 are identified 
as sunspots. The remaining active pixels are taken as network 
and faculae, with the following exclusions (see detailed discus- 
sion in Yeal3): 


— Yeal3 noted that the magnetic canopy of sunspots extend 
beyond their boundary in the continuum image. As in that 
study, we excluded the active pixels contiguous to sunspots. 
Although some of these points might correspond to network 
and faculae, the apparent magnetic flux density is affected 
by the encroaching magnetic canopy and stray light originat- 
ing from the sunspots. This measure also accounts for bright 
penumbra filaments that are not captured by the continuum 
intensity criterion. 

— Stand-alone active pixels are excluded to minimise the incor- 
rect inclusion of magnetogram noise as network and faculae. 

— Active pixels where B > 800 G are excluded because 
they mainly correspond to bright features associated with 
sunspots near the limb than network and faculae. See the de- 
tailed discussion in Yeal3. 

— Close to the limb, the B = Bj,/u approximation breaks down 
from the combination of magnetogram noise and foreshort- 
ening. For this reason, we ignored the active pixels where 
u < 0.1. The effect on the study is very light because only 
1% of the solar disc, by area, lies outside the u = 0.1 loci. 


It is worth emphasising here that it is not necessary in such an 
investigation to include the network and faculae as completely 
as possible but rather to minimise false positives. While these 
measures inevitably result in some network and faculae being 
excluded, they minimise false positives while still retaining most 
of the network and faculae (Fig. 1). A total of 5.0 x 107 image 
pixels are identified as corresponding to network and faculae. 
They span By, < B < 800 G and 0.1 < u < 1.0, the widest range 
of magnetic flux densities and disc positions we can reasonably 
examine with HMI observations. 
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c) 


Fig. 1. Top: HMI and AIA observations from one of the data days, 12 June 2014. From left to right, the intensity image at a) the 
continuum and b) the core of the Fe I 6173 À line, and c) the AIA 1700 À channel, d) the longitudinal magnetogram, and e) the 
map of the image pixels identified as network and faculae (white). The image pixels identified as sunspots are marked in red. The 
magnetic activity contiguous to sunspots and pores, marked in blue, is excluded from the analysis (see Sect. 3.1). The greyscale is 
saturated for the intensity images at null and 150% of the quiet-Sun level at disc centre, and at +50 G for the magnetogram. Bottom: 
Corresponding zoomed-in inset of the boxed 300 x 300 arcsec region. 


3.2. Intensity contrast 


Following Yeal3, the intensity contrast at position (x, y) in a 
given image C (x, y) is defined as 


ege 


= -1 1 
Tgsn (x, y) ' t ) 


where J and /4,, denote the apparent intensity and the quiet-Sun 
level, respectively. 

Taking each intensity image, we masked the points where 
B > 10 G in the corresponding longitudinal magnetogram, leav- 
ing the quiet Sun. The quiet-Sun intensity as a function of dis- 
tance from disc centre, that is, the limb-darkening profile, is 
given by the fifth-order polynomial in u fit to this quiet-Sun 
image (following Neckel & Labs 1994). We normalised the 
quiet-Sun image by the limb-darkening profile. Image distor- 
tions from inhomogeneities in instrumental properties across the 
image plane manifest in the normalised quiet-Sun image as devi- 
ations in the overall level from unity. The flat-field image, map- 
ping these image distortions, is given by the polynomial surface 
fit to the normalised quiet-Sun image (after first smoothing it 
with a 401 x 401 pixel boxcar filter to filter out small-scale vari- 
ability). For the 1700 A images, where the image distortions are 
most severe, the root mean square (RMS) difference between 
the flat-field image and unity ranged from 3.2% to 4.6% over 
the data days. We corrected the original intensity image for im- 
age distortions by dividing it by the flat-field image before cal- 
culating the intensity contrast (Equation 1). We determined the 
limb-darkening profile and flat-field image of each intensity im- 
age separately on the observation that HMI and AJA instrument 
throughput and the image distortions appear to vary with time. 

The line core intensity is modulated not just by fluctuations 
in line strength and shape, but also by continuum excess that fills 
in the spectral line. To aid the interpretation of line core intensity 
contrasts, we factored out continuum changes by dividing each 
line core intensity image by the corresponding normalised con- 


tinuum intensity image before we calculated the intensity con- 
trast. (With this normalisation, the line core intensity image in- 
dicates the level if the continuum is fixed at the quiet-Sun level.) 
We refer to Yeal3, who examined the effect of this step on ap- 
parent line core intensity contrasts in detail. 


The data reduction (Sect. 2.3) and the process by which we 
isolated network and faculae (Sect. 3.1), and determined their in- 
tensity contrast, just described, closely follows Yeal3. The only 
significant departure from Yeal3 is that we resampled the HMI 
observations to the image pixel scale of AIA in order to extend 
the study to include the AIA 1700 À channel (Sect. 2.3). As we 
demonstrate in Sect. 4.1.1, this had negligible effect on the ap- 
parent intensity contrast in the HMI passbands. 


3.8. Local magnetic flux density 


One of the aims of this study is to examine the effect of how 
magnetically crowded a given region of the solar surface is on 
the radiant properties of the enclosed network and faculae. To 
this end, we studied the variation in intensity contrast with what 
we call the local magnetic flux density, (B)j4,, which is defined 
as the mean absolute magnetic flux density in the 9 x 9 pixel 
(5.4 x 5.4 arcsec) window centred on each network and facular 
image pixel. To minimise any bias from magnetogram noise, we 
smoothed each longitudinal magnetogram with a Gaussian ker- 
nel with an FWHM of 3 pixels before we calculated this quantity. 


We adopted a window size of 9 x 9 pixel. Tests showed that 
the window size has no qualitative effect on the results, to be 
presented in Sect. 4.3. The effect of magnetic concentrations on 
one another evidently diminishes with separation such that only 
the magnetic fields closest to a given point is relevant. Therefore 
broadening the window beyond a certain minimum does not add 
more of the relevant neighbouring magnetic activity, but merely 
scatters the resulting (B);.. 
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4. Results 
4.1. Variation with disc position and magnetic flux density 


First, we examine the variation in intensity contrast with distance 
from disc centre and magnetic flux density, represented by u and 
B, respectively. The intensity contrast measurements cover 0.1 « 
H € 1.0and Bnr < B < 800 G. To elucidate the CLV, we grouped 
the measurements by B into eight intervals and charted the trend 
with u within each interval (Fig. 2). As for the B dependence, we 
divided the solar disc into eight concentric annuli of equal area 
and examined the variation in intensity contrast with B within 
each annulus (Fig. 3). 


4.1.1. Empirical relationship 


Following Ortiz et al. (2002) and Yeal3, we also encapsulated 
the CLV and B dependence of intensity contrast at each passband 
into an empirical relationship given by the bivariate polynomial 
in u and B fit to the intensity contrast measurements. We binned 
the intensity contrast measurements by u and In B in intervals of 
0.02 and 0.05, respectively. We then fit a cubic polynomial in u 
and B to the bin-averaged intensity contrast. In other words, we 
described the intensity contrast by a relationship of the form 


0 47 0 
B do) do do aos || H 
4i0 411 412 413 H 
C(u, B) = B 2] (2) 
420 421 022 423 H 
B? 430 431 432 433 we 


where aj; denotes the coefficient of the Biu’ term, tabulated in 
Table 1. This captures the variation in intensity contrast with 4 
and B, as demonstrated by the alignment to the measured inten- 
sity contrast (Figs. 2 and 3). 

In order to include the AIA 1700 À channel in this study, we 
resampled the HMI observations to the image pixel scale of AIA 
(Sect. 2.3). Figs. 2a and 3a show that the fit to the measured in- 
tensity contrast at the continuum of the Fe I 6173 A line (solid 
curves) is almost exactly similar to that from the same analysis 
in Yeal3 with HMI observations at the original spatial resolution 
(dashed curves). We noted the same at the Fe I 6173 A line core, 
which is not shown to avoid cluttering the plots. Resampling the 
HMI observations to the AIA image pixel scale had no apprecia- 
ble effect on the apparent intensity contrasts. 


4.1.2. Fe | 6173 A line continuum 


The continuum of the Fe I 6173 A line is formed in the lower 
photosphere, where magnetic flux tubes are heated through the 
side walls of the opacity depression by radiation from the sur- 
rounding convection (Spruit 1976). The intensity of network and 
faculae is modulated by the competing influence of this lateral 
heating and the magnetic suppression of convection within mag- 
netic concentrations (Spruit & Zwaan 1981; Grossmann-Doerth 
et al. 1994). 

With increasing distance from disc centre to limb, the inten- 
sity contrast increases to a maximum before declining towards 
the limb (Fig. 2a). The initial climb corresponds to the hot side 
wall of the opacity depression rotating into view and the cool 
centre becoming increasingly occulted. The side wall is eventu- 
ally blocked as well, giving the drop near the limb. The aver- 
age magnetic flux tube size increases with B (Ortiz et al. 2002; 
Yea13), with the consequence that lateral heating becomes less 
and less efficient. This contributed to the observation here that 


the slope of the contrast-versus-B profile diminishes with B (Fig. 
3a). 


4.1.3. Fe 16173 A line core and 1700 A 


The core of the Fe I 6173 A line and the AIA 1700 A chan- 
nel are formed over the middle photosphere, where magnetic 
flux tubes are heated by radiation from deeper layers (Knólker 
et al. 1991), and mechanical and Ohmic dissipations (Moll et al. 
2012). Mechanical heating includes longitudinal tube waves 
(LTWs) that are excited in magnetic flux tubes by surround- 
ing convection (Carlsson & Stein 1992; Fawzy et al. 2012). The 
LTWSs steepen as they propagate upwards, eventually forming 
shocks and dissipating their energy. This is one avenue through 
which the temperature structure and therefore the radiant proper- 
ties of magnetic flux tubes can be influenced by the inhibition of 
convection in their surroundings by nearby magnetic fields (c.f. 
Sect. 4.3). 

The Fe I 6173 À line core intensity is also modulated by 
Zeeman line changes and continuum excess filling in the spectra 
line, although the latter is factored out here in the derivation of 
intensity contrast (Sect. 3.2). The effect of Zeeman line changes 
on intensity contrast is indicated by the following: The HMI line 
width maps concurrent to the HMI observations examined here 
are available for 27 of the 100 data days. Taking these data, we 
found that the line width at network and faculae increases mono- 
tonically with B, depicted for disc centre (0.94 < u < 1.00) in 
Fig. 4. Because the Zeeman effect does not change the equiva- 
lent width of a given spectral line (except if itis saturated, which 
is not the case here), an increase in line width has to be accom- 
panied by a decrease in line depth. Therefore, the steady rise in 
line width with B indicates a similar trend in line core intensity. 

As opposed to what we observed at the continuum of the 
Fe I 6173 À line (Fig. 2a), at these two passbands, the inten- 
sity contrast declines steadily from disc centre to limb (Fig. 2b). 
We quantify the strength of this decline by the proportional drop 
between disc centre and u = 0.2, 1 — C (u = 0.2) /C (u = 1.0) 
(Fig. 5). For each interval of B in Fig. 2b, we determined 
1 - € (u = 0.2) /C (u = 1.0) from the cubic polynomial fit to the 
CLV profile (except for the By, < B < 50 G interval, as no data 
are available around u = 0.2 here). The uncertainty is propagated 
from the regression error. We also estimated this quantity from 
the empirical intensity contrast models (Equation 2 and Table 1). 

At 1700 A (blue, Fig. 5), 1-C (u = 0.2) /C (u = 1.0) declines 
steadily with B, while at the Fe I 6173 A line core (red), this 
quantity falls less steeply with B and only up to about 400 G 
before it rises instead. The following physical processes are at 
work: 


1. Magnetic flux tubes tend towards a surface normal orienta- 
tion. This means, going from disc centre to the limb, rays 
change from passing through single flux tubes to passing in 
and out of multiple flux tubes (Biinte et al. 1993; Solanki 
et al. 1998). This leads to a drop in intensity contrast as radi- 
ation from flux tubes becomes absorbed into the interleaving 
non-magnetised atmosphere that is traversed by the rays. Of 
course, the greater the magnetic filling factor, the weaker this 
effect. 

2. The temperature gradient of the solar atmosphere is weaker 
within magnetic concentrations (e.g. Vernazza et al. 1981; 
Fontenla et al. 1991) such that the temperature contrast in- 
creases with height. As both the Fe 16173 A line and the AIA 
1700 A channel are optically thick, going from disc centre to 
limb, we look shallower and shallower into the solar atmo- 
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Fig. 2. a) Intensity contrast at the continuum of the Fe I 6173 A line in the annotated intervals of B. We binned the intensity contrast 
measurements within each B interval by u in bins of 0.05. The filled circles and error bars correspond to the mean and standard 
deviation within each u bin. The standard error of the mean is given by the quotient of the standard deviation and the square 
root of the number of measurements in the bin. Each bin contains 10? measurements, such that the standard error is too small to 
be represented in this plot. The solid curves follow the bivariate polynomial in the u and B fit to the measured intensity contrast 
(Equation 2 and Table 1). The dashed curves represent the result from the similar analysis in Yeal3 at the original spatial resolution 
of HMI, which are largely hidden because they agree so closely. b) The same, except for the intensity contrast at the core of the Fe 
1 6173 A line (red) and 1700 A (blue). The Fe I 6173 A line core values are scaled by a factor of 10 to facilitate the comparison. 
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Fig. 3. a) Intensity contrast at the continuum of the Fe I 6173 A line as a function of B from about disc centre (0.94 < u < 1.00) to 
near the limb (0.10 < u < 0.36). At each interval of u, we binned the intensity contrast measurements by B in 40 G bins. The mean 
and standard deviation within each B bin is represented by the filled circles and error bars. As in Fig. 2, the standard error of the 
mean is negligible because of the size of our data set. The solid curves represent the bivariate polynomial in u and B fit to measured 
intensity contrast (Equation 2 and Table 1) and the dashed curves show the same from Yeal3, which is almost indistinguishable 
because of the close agreement. b) The corresponding plots for the core of the Fe I 6173 A line (red) and 1700 A (blue). Again 
following Fig. 2, we scaled the Fe I 6173 A line core values by a factor of 10. 
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Table 1. Coefficients of the bivariate polynomial in the u and B (Equation 2) fit to the measured intensity contrast. 
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Fig. 4. As a function of B, the width of the Fe I 6173 Å line, 
relative to the quiet-Sun level, about disc centre (0.94 < u < 
1.00). The relative line width of network and faculae, which is 
essentially their contrast in the line width map, is derived in an 
analogous manner to the intensity contrast (Sect. 3.2). The filled 
circles and error bars have the same meaning as in Fig. 3. 
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Fig. 5. Proportional drop in intensity contrast between disc cen- 
tre and u = 02, 1 — C (u = 0.2) /C (u = 1.0) at the Fe I 6173 
À line core (red) and 1700 A (blue). The plot points denote the 
estimate from the intensity contrast CLV profiles in Fig. 2b. The 
curves follow the solution from the empirical intensity contrast 
models (Equation 2 and Table 1). The dashed line marks the min- 
imum of the Fe I 6173 A line core curve. 


sphere, resulting in an increase in temperature and intensity 
contrast. Obviously, this effect becomes stronger with mag- 
netic filling factor. This mechanism counters the decline in 
intensity contrast from disc centre to limb (Fig. 2b) such that 
it is less steep than it would be without this effect. 

. The Fe I 6173 À line weakens from disc centre to the limb 
(Norton et al. 2006), which also decreases the contribution 
by Zeeman line changes to the intensity contrast. As Zeeman 
line changes become stronger with B, so does the disc centre 
to limb decline in intensity contrast from this effect. 


As B and therefore the magnetic filling factor rises, the first two 
mechanisms lead to a weaker disc centre-to-limb decline in in- 
tensity contrast, while the third mechanism, which is only ap- 
plicable to the Fe I 6173 À line core, does the opposite. The 
action of the first two mechanisms at 1700 À produces the 
monotonic decline in 1 — C (u = 0.2) /C (u = 1.0) with B. The 
competition between the first two mechanisms and the third 
at the Fe I 6173 À line core produces the reversal in trend in 
1 - C (u = 0.2) /C (u = 1.0) with B at around 400 G. 

At about the disc centre (0.94 < u < 1.00, Fig. 3b), the 
intensity contrast rises with B but at a diminishing rate. While 
the intensity contrast at the Fe I 6173 À line core (red) increases 
monotonically with B, at 1700 À (blue) it peaks at about 350 
G before it starts to decline. We explain the diminishing slope 
and the divergent trend above 350 G next in Sect. 4.2. From disc 
centre (0.94 < u < 1.00, Fig. 3b) to the limb (0.10 < u < 
0.36), the B dependence of the intensity contrast changes. At 
the Fe I 6173 A line core (red), the intensity contrast changes 
from increasing monotonically with B to peaking at intermediate 
levels of B, while the converse is apparent at 1700 A (blue). The 
reason is the B dependence of the CLV at the two passbands (Fig. 
5), as discussed above. 


4.2. Heating efficiency of the magnetic flux 


We denote with Umax(B) and Cmax(B) the position of and inten- 
sity contrast at the peak of the intensity contrast CLV profile at 
a given level of B. Following Ortiz et al. (2002) and Yeal3, we 
examined the heating efficiency of magnetic flux as indicated by 
the intensity contrast per unit magnetic flux density at the CLV 
profile peak, Cmax/B, which we call the specific contrast. 

At each passband and each interval of B in Fig. 2, we esti- 
mated Umax and Cmax from the cubic polynomial fit to the CLV 
profile. The specific contrast is then given by the quotient of C max 
and the mean B of the network and facular points in the B in- 
terval. The corresponding uncertainty is propagated from the re- 
gression error and the standard error of the mean B. We excluded 
the Bg, < B € 50 G interval, where this analysis is adversely af- 
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Fig. 6. a) As a function of B, specific contrast, Cmax/B (black, 
left-hand axis) and the position of Cmax, Umax (grey, right-hand 
axis) at the continuum of the Fe I 6173 À line. The plot points 
correspond to the values derived from the CLV profiles (Fig. 2a) 
and the curves to the estimate from the empirical intensity con- 
trast models (Equation 2 and Table 1). b) Specific contrast at the 
Fe I 6173 À line core (red) and 1700 A (blue). To facilitate the 
comparison, we scaled the Fe I 6173 A line core values by a 
factor of 11.9 to match the 1700 A values at the first plot point. 
At these two passbands, Umax = 1 at all levels of B, hence the 
omission from the figure. 


fected by the fact that the apparent intensity contrast and B are 
weakest and therefore most uncertain here, and extend a limited 
range of u. We also determined max and Cmax/B from the em- 
pirical intensity contrast models (Equation 2 and Table 1). The 
results are depicted in Fig. 6. 

At the continuum of the Fe I 6173 À line, the intensity con- 
trast peaks closer and closer to the limb with increasing B (grey, 
Fig. 6a). This was similarly noted by Ortiz et al. (2002) and 
Yea13, who asserted, based on the hot wall model (Spruit 1976), 
that this indicates the size of magnetic flux tubes on average in- 
creases with B. At the two other passbands, {max = 1.0 at all lev- 
els of B (Fig. 2b). The divergence stems from the difference in 
formation height and therefore the physical processes that drive 
the magnetic enhancement of intensity, as discussed in Sect. 4.1. 

The specific contrast at the continuum of the Fe 16173 A line 
(black, Fig. 6a) rises with B, up to a maximum at about 150 G, 
before it starts to decline. At the Fe I 6173 A line core (red, Fig. 
6b) and 1700 À (blue), the specific contrast, given here by the 
disc centre intensity contrast per unit B, declines monotonically 
with B and more rapidly so at the latter. As noted in Sect. 4.1.3, 
around disc centre (0.94 < u x 1.00, Fig. 3b), the Fe I 6173 A 
line core (red) and the 1700 A intensity contrast (blue) therefore 
rises with B at a diminishing rate and more rapidly so at 1700 A, 
such that the trend with B reverses at around 350 G. 


As argued by Yeal3, in non-sunspot magnetic features in 
full-disc observations, the points with lower magnetogram sig- 
nals mainly correspond to the network, and the higher levels 
correspond to faculae. Therefore, the decline in specific contrast 
with B indicates that the quiet-Sun network for the same amount 
of magnetic flux is brighter and hotter than active region fac- 
ulae. For a given magnetic flux tube, radiative heating through 
the side walls of the opacity depression and mechanical heating 
by LTWSs are inhibited when convection in the greater area is 
suppressed by nearby magnetic fields. Because facular regions 
evidently contain more magnetic fields than the quiet Sun, this 
effect is greater on faculae than in the network. We examine this 
in greater detail in the following section (Sect. 4.3). Another fac- 
tor that could have contributed to the downward trend in specific 
contrast with B at the continuum of the Fe I 6173 À line is the 
effect of the increase in average flux tube size with B on the effi- 
ciency of the lateral radiative heating. 

The divergent B dependence at the continuum and core of 
the Fe I 6173 A line below 150 G was similarly noted by Yeal3, 
who attributed it to the tendency of magnetic concentrations to 
aggregate in dark intergranular lanes on the intensity contrast in 
the lower photosphere, that is, the continuum (Title & Berger 
1996; Schnerr & Spruit 2011). The steeper decline with B at 
1700 À as compared to the Fe I 6173 À line core is at least 
partly due to the contribution by the Zeeman line changes to 
the intensity contrast at the latter, which increases with B (Sect. 
4.1.3). 


4.3. Effect of the local magnetic flux 


In this section, we examine the effect on the intensity contrast of 
the network and faculae through nearby magnetic fields, which 
are represented here by the local magnetic flux density, (Boc 
(defined in Sect. 3.3). To this end, we re-examine the CLV at 
low and high B (Fig. 7a), the B dependence near disc centre and 
the limb (Fig. 7b), and the specific contrast (Fig. 8) at the points 
where (Boc € 100 G (solid curves) and (Boc > 100 G (dotted 
curves) separately. We then compare the results to what we de- 
rived earlier in Sects. 4.1 and 4.2, for which we take all available 
data together (filled circles). At lower B, the values from the en- 
semble lie close to the (Boc < 100 G values and at higher B to 
the (B), > 100 G values. This reflects the fact that most of the 
points with lower (higher) B are in the less (more) magnetically 
crowded network (faculae) fields. 

First, we consider the continuum of the Fe I 6173 À line 
(black, Figs. 7 and 8). At low B (50 G < B < 100 G, Fig. 
7a), the intensity contrast increases from (B), < 100 G to 
(B}ioc > 100 G and peaks closer to the limb. At high B (500 G < 
B < 640 G, Fig. 7a), the intensity contrast is suppressed. The 
CLV profile peak is at about the same position, but the decline 
limb-wards of the peak is weaker such that the intensity contrast 
is higher near the limb. The effect of nearby magnetic fields on 
the intensity contrast varies gradually with B between what is 
noted here at the lower and upper bounds, not shown for brevity. 
When we consider that the CLV is from the heated side walls 
of magnetic flux tubes coming into greater view before being 
blocked by the side facing away from the observer (cf. Sect. 
4.1.2), the limb-ward shift of the CLV profile peak at low B and 
the intensity contrast enhancement near the limb at high B indi- 
cate greater flux tube sizes in magnetically crowded regions. 

Most of the studies of the continuum intensity contrast of 
the network and faculae at disc centre found that it declines into 
the negative as B approaches 0 G (see Kahil et al. 2016, and the 
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Fig. 7. a) CLV of the intensity contrast at low (left) and high B (right), at (B), < 100 G (solid curves) and (B), > 100 G (dotted 
curves). From top to bottom, we show the values from the continuum (black) and core (red) of the Fe I 6173 Å line and 1700 Å 
(blue). b) Variation in intensity contrast with B about disc centre (left) and near the limb (right), in the same (B),,, ranges. The 
intensity contrast profiles are derived in the same manner as the corresponding profiles, which were computed taking all available 
data into account (filled circles, taken from Figs. 2 and 3). As in Figs. 2 and 3, the standard error is negligible because of the size of 
the data set. 
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Fig. 8. Same as Fig. 6, except at (B), x 100 G (solid curves) 
and (Boc > 100 G (dotted curves). The filled circles denote the 
values from taking all available data into consideration, adopted 
from the earlier figure. The uncertainty of the (B), < 100 G and 
(B). > 100 G series is of similar order as the ensemble series. 
It is omitted here to aid visibility. 


summary of earlier results by Yea13). Both here (Fig. 6) and in 
Yea13, we find that below a certain level of B, the specific con- 
trast in the continuum starts to decline. Both observations have 
been attributed to the tendency for magnetic concentrations to 
aggregate in the dark intergranular lanes (see the detailed ar- 
gument in Yeal3). Fig. 7b shows that as B approaches 0 G, 
the intensity contrast is enhanced by the nearby magnetic fields. 
Below 150 G, the specific contrast at (B), < 100 G (solid curve, 
Fig. 8a) declines, as seen in the ensemble (filled circles), but the 
specific contrast at (Boc > 100 G (dotted curve) rises instead. 
We surmise that in magnetically crowded regions, the magnetic 
concentrations are overall less deeply embedded in intergranular 
lanes than in quieter regions. 


Next we examine the Fe I 6173 À line core (red, Figs. 7 and 
8) and 1700 A (blue). At all levels of B, the intensity contrast 
from disc centre to limb changes from being suppressed to being 
enhanced by the local magnetic flux (Fig. 7b), such that the CLV 
is weaker (Fig. 7a). As noted in Sect. 4.1.3, the decline in inten- 
sity contrast from disc centre to the limb is in part due to inten- 
sity excess produced in magnetic flux tubes becomes absorbed 
into the interleaving non-magnetised atmosphere between flux 
tubes. For a given magnetic filling factor, the larger the flux 
tubes, the weaker this effect (Solanki et al. 1998). The weaker 
disc centre-to-limb drop is therefore another indication of greater 
flux tube sizes in magnetically crowded regions. In combination 
with the lower heating efficiency (Fig. 8b), the intensity contrast 
switches from being suppressed to being enhanced by the local 
magnetic flux from disc centre to limb. 


At each passband, except below about 100 G in the contin- 
uum of the Fe I 6173 A line, the specific contrast is higher at 
(B). < 100 G (solid curves, Fig. 8) than at (By, > 100 G 
(dotted curves). As noted at the beginning of this section, from 
low to high B, the values calculated with the ensemble (filled 
circles) migrate from the level at (B), < 100 G to the level 
at (B), > 100 G because most of the lower (higher) B points 
came from magnetically sparsely (crowded) regions. As asserted 
in Sect. 4.2, the downward trend in heating efficiency with B 
seen in the ensemble is partly due to the growing prevalence and 
influence of nearby magnetic fields. 

The analysis here revealed that even though the heating effi- 
ciency is largely suppressed by nearby magnetic fields, intensity 
contrast can be enhanced instead at certain B and yw. This is due 
to the effect of how magnetically crowded a given region is on 
flux tube sizes and how deeply embedded magnetic concentra- 
tions are in intergranular lanes. 


4.4. Comparison of the quiet-Sun network and active region 
faculae 


Next, we investigate if the differences in intensity contrast be- 
tween the quiet-Sun network and active region faculae noted by 
Lawrence et al. (1993) and Kobel et al. (2011), discussed in the 
introduction, can be seen in the present data set, and if it can be 
accounted for by the effect of nearby magnetic fields on the in- 
tensity contrast alone. The implication being that if this were not 
the case, it implies that there is a fundamental physical difference 
between the quiet-Sun network and active region faculae. 

For each data day, we segmented the solar disc into the 
quiet Sun and active regions, broadly following Hagenaar et al. 
(2003). We refer to any contiguous patch of magnetically active 
pixels (i.e. where Bios > 308) as a magnetic island. We isolated 
the magnetic islands where Bio, > 50g, over an area of at least 
11x11 pixels (6.6x 6.6 arcsec). Taking these magnetic islands to 
be the core of active regions, we expanded their boundary with a 
51x51 pixel (30.6 x 30.6 arcsec) boxcar filter. Taking the points 
isolated in Sect. 3.1, we counted the points within the boundary 
as active region faculae and outside the boundary as quiet-Sun 
network. An example of the segmentation is depicted in Fig. 9a. 
As expected, the distribution of (B); is higher in active region 
faculae than in the quiet-Sun network (Fig. 9b). 

Following Kobel et al. (2011), we compared the continuum 
intensity contrast in the quiet Sun (purple, Figs. 9c and 9d) and 
in active regions (green) at disc centre (0.99 < u < 1.00). The 
result is qualitatively similar to what was reported by the earlier 
study; compared to active region faculae, the quiet-Sun network 
intensity contrast is higher about the peak of the contrast versus 
B profile, and lower approaching B = 0 G. Taking all the points 
identified as network and faculae at 0.99 < u < 1.00, we binned 
the points by B and (B), (in intervals of 20 G) and fit a quartic 
bivariate polynomial in B and (B),,, to the bin-averaged intensity 
contrast. The fit is in effect an empirical model of disc-centre in- 
tensity contrast as a function of B and (B),,, alone (i.e. assuming 
no fundamental physical difference between network and facu- 
lae). We re-derived the contrast versus B profiles by calculating 
the intensity contrast of each network and facular image pixel 
from its B and (Boc with this empirical model instead of taking 
the apparent contrast in the intensity image. As evident in Figs. 
9c and 9d, the modelled profiles (curves) are in close agreement 
with the measured profiles (filled circles). In Fig. 10 we show 
similar results from repeating this analysis on the Fe I 6173 A 
line core and 1700 A. The results here suggest that the apparent 
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Fig. 9. a) Taken from Fig. 1, the HMI longitudi- 
nal magnetogram from 12 June 2014. The yel- 
low contours encompass the areas we identified 
as active regions. For the entire data set (i.e. all 
data days), b) (Boc and c) the Fe 16173 À con- 
tinuum intensity contrast of the quiet-Sun net- 
work (purple) and active region faculae (green) 
at disc centre (0.99 < u < 1.00). The filled cir- 
cles represent the mean within successive in- 
tervals of 40 G. In b) the error bars denote 
the corresponding standard deviation, and in c) 
they denote the standard error. The curves cor- 
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Fig. 10. Similar to Fig. 9c, except for a) the Fe 
16173 A line core and b) 1700 A. As in the ear- 
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divergence between the quiet-Sun network and active region fac- 
ulae intensity contrast arises is because active regions are more 
magnetically crowded (Fig. 9b). If the apparent divergence be- 
tween quiet Sun and active region intensity contrast were statis- 
tical or would arise from factors other than the local magnetic 
flux density, then the modelled profiles should overlap with one 
another instead of reproducing the difference between the mea- 
sured quiet-Sun and active region profiles. This is supported by 
the following. 


We binned the points by B, (Boc , and u (in intervals of 
20 G, 20 G, and 0.25, respectively) separately for the quiet-Sun 
network and active region faculae, and took the bin-averaged in- 
tensity contrast. The scatter plot of active region faculae versus 
quiet-Sun network bin-averaged intensity contrast at each pass- 
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lier figure, the filled circles and error bars rep- 
resent the mean and standard error within suc- 
cessive intervals of 40 G. 


band is depicted in Fig. 11. The correlation (Pearson’s R?) be- 
tween network and facular intensity contrast and the slope of 
the straight-line fit to the scatter plot, indicated in the figure, are 
within about 1% to 4% of unity. In other words, after factoring 
out differences in B, (Boc , and p, the network and facular in- 
tensity contrast are similar to this margin. 


There is no indication, at least within the limits of this study, 
of any fundamental physical difference between network and 
faculae. The different radiant behaviour of the quiet-Sun net- 
work and active region faculae arise because the faculae face 
more nearby magnetic fields alone. The results here and in Sect. 
4.3 highlight the need for models of solar irradiance to include 
the effect of nearby magnetic fields on the network and facular 
intensity. This is so far lacking in existing models. 
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Fig. 11. Scatter plot of active region faculae vs. quiet-Sun network intensity contrast at a) the continuum and b) the core of the Fe I 
6173 Å line, and c) at 1700 Å. The plot points represent the result of binning the intensity contrast measurements by B, (Boc , and 
4, and taking the average within each bin (see Sect. 4.4). The straight lines denote the linear fit to each scatter plot, constrained to 
pass through the origin. The Pearson correlation coefficient R? and the slope of the linear fit are annotated. 


4.5. Comparison of the Fe | 6173 A line core and 1700 A 


Here, we investigate what drives the divergence between in- 
tensity contrast at the Fe I 6173 À line core and 1700 À in 
addition to the Zeeman line changes in the former (discussed 
in Sects. 4.1 and 4.2). We study the intensity contrast CLV at 
50 G < B < 100 G (Fig. 12a). Assuming the solar atmosphere 
radiates like a black body, we estimated from the measured in- 
tensity contrast at 1700 À (blue circles) the intensity contrast at 
the Fe I 6173 À line core (blue curve). We also projected the 
measured intensity contrast at the Fe I 6173 À line core (red cir- 
cles) to 1700 A (red curve). This calculation requires as input the 
effective temperature of the quiet Sun at both passbands. We as- 
sumed an effective temperature of 5800 K at the Fe I 6173 A line 
core. We then set the effective temperature at 1700 À at 6020 K; 
this is the level that optimises the agreement between measured 
and projected intensity contrast. (The effective temperature of 
the quiet Sun at the two passbands cannot be identical given the 
difference in formation height, cf. Sect. 2.) This analysis, how- 
ever approximate, shows an excellent agreement between mea- 
sured and projected intensity contrast. 

We confine this discussion to the low B regime because this 
is where the Zeeman line changes are weakest and can be fac- 
tored out. For information, we depict the results of repeating the 
above calculation at high B in Fig. 2b. As expected, the pro- 
jection, which ignores the Zeeman line changes, fails to repro- 
duce the measured intensity contrast. In addition to the Zeeman 
line changes, the divergence between the intensity contrast at 
the Fe I 6173 À line core and 1700 À arises from the different 
wavelength and temperature regime (from the different forma- 
tion height), as governed by the black-body law. 

Photospheric iron lines, including the Fe I 6173 Á, are 
formed largely in LTE. Departures from LTE are sufficiently 
minute such that observed line profiles can be equally well repli- 
cated by LTE and non-LTE radiative transfer schemes with an 
adjustment of the iron abundance (Rutten & Kostik 1982; Bruls 
et al. 1991; Shchukina & Trujillo Bueno 2001). On the other 
hand, it is not understood whether the AIA 1700 À channel is 
formed in LTE or non-LTE, due in part to the forbiddingly com- 


plicated zoo of relevant species (see, e.g. Fossum & Carlsson 
2005). Therefore, the fact that we can replicate intensity contrast 
at one passband with the measurements from the other assuming 
black-body behaviour suggests either that the AIA 1700 A chan- 
nel is formed in LTE or that the AIA 1700 À channel is formed 
in non-LTE, but departures from LTE affect the non-magnetised 
and magnetised solar atmosphere in such a way that the intensity 
contrast is not affected. As the effective temperature at 1700 Ais 
a fit parameter, set at the level that optimises the agreement be- 
tween measured and modelled intensity contrast, the agreement 
between the two could also mean that the AIA 1700 A channel 
is formed in non-LTE and non-LTE effects on the intensity con- 
trast can be reproduced by assuming the appropriate effective 
temperature. 

If indeed departures from LTE affect the non-magnetised and 
magnetised solar atmosphere in such a way that the intensity 
contrast is unaffected, it will explain the following: Certain mod- 
els of solar irradiance (Krivova et al. 2006; Ball et al. 2011; Yeo 
et al. 2014) make use of the intensity contrast of solar surface 
magnetic features that are calculated from model solar atmo- 
spheres assuming LTE (Unruh et al. 1999). The Yeo et al. (2014) 
model reproduces the observed variability in the wavelength- 
resolved or spectral solar irradiance all the way down to about 
160 nm even though the absolute level already starts to differ 
from measurements at below 300 nm as the LTE assumption 
breaks down (Yeo et al. (2014). 


4.6. Comparison with the Ca II K line 


Various studies have reported power-law relationships between 
chromospheric emission and magnetic flux density (Schrijver 
et al. 1989; Harvey & White 1999; Ortiz & Rast 2005; Rezaei 
et al. 2007; Loukitcheva et al. 2009; Kahil et al. 2016). Notably, 
the study by Harvey & White (1999), the only to cover the full 
solar disc, found the relationship between intensity excess at the 
Ca II K line, the difference in intensity to the quiet-Sun level, 
and Bios to be invariant with u. We investigate if emission at the 
passbands examined in this study, formed in the lower and mid- 
dle photosphere (Sect. 2), conform to similar relationships with 
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Fig. 12. Filled circles: Intensity contrast CLV at the Fe I 6173 A line core (red) and 1700 A (blue), at a) 50 G < B < 100 Gand b) 
640 G < B < 800 G, taken from Fig. 2b. Curves: Result of estimating the intensity contrast at one passband from the measurements 
at the other, assuming the solar atmosphere radiates like a black body (solid) and the corresponding 1c bounds from the uncertainty 


in the effective temperature at 1700 A (dashed). See Sect. 4.5. 
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Fig. 13. Log-log scatter plot of the intensity excess and Bios at a) the continuum and b) the core of the Fe I 6173 À line, and c) 1700 
A. The intensity excess is given in units of ppm of the disc-integrated quiet-Sun intensity. The plot points, colour-coded by u, denote 
the result of grouping the network and facular image pixels with similar Bios and u together, and taking the average intensity excess. 


See Sect. 4.6. 


Bios and u. In keeping with the cited studies, we study the rela- 
tionship between intensity excess, denoted AJ, and Bios, instead 
of intensity contrast and B, as in the rest of this study. 


In Fig. 13 we show the log-log scatter plot of AJ and Bios 
at each passband. A power-law relationship between AJ and Bios 
that is invariant with u will show itself in such a representation as 
the plot points, coming from across the entire disc, falling along 
the same straight line. As is evident from the figure, AJ does not 
conform to a power-law relationship with Bios, and neither is the 
relationship with Bios invariant with yz at any of the passbands. 


What is the implication for models of solar irradiance that 
use chromospheric indices, including indices based on the Ca II 
K line (Keil et al. 1998), as direct indications of the effect of net- 
work and facular on solar irradiance? The use of chromospheric 
indices as proxies of network and facular brightening is based 
on the fact that they are dominantly modulated by the chromo- 
spheric network and plage overlaying photospheric network and 
faculae. With this in consideration, we examine the variation in 
disc-integrated intensity, from the changing prevalence and dis- 
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tribution of network and faculae, at the passbands examined here 
and at the Ca II K line. 


The HMI has been in operation since 30 April 2010. For 
each day since, up to 31 October 2018, we took a 315s longi- 
tudinal magnetogram and continuum intensity image, and iden- 
tified the network and faculae as done in Fig. 3.1, except that we 
did not exclude the non-sunspot magnetic activity contiguous to 
sunspots and pores (this makes no difference to the current anal- 
ysis). At each passband, using the empirical model of intensity 
contrast as a function of B and u derived in Sect. 4.1.1, we de- 
termined the daily disc-integrated network and facular intensity 
excess (black, Fig. 14a). To elucidate the variability over solar 
cycle and rotation timescales, we also draw the 81-day moving 
average (Fig. 14b) and the zoomed inset of 2014 (Fig. 14c). 


Harvey & White (1999) reported a power-law exponent of 
0.5, or in other words, AJ œ «Bi, at the Ca II K line. For each 
day, we took the sum of Bi, of the network and facular image 
pixels. The result is a time series with the variability, although 
not the absolute scale, of the disc-integrated network and facu- 
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Fig. 14. a) Black: Daily disc-integrated intensity excess at the continuum and core of the Fe I 6173 A line and 1700 A (top to 
bottom), calculated with the empirical intensity contrast model derived in this study (Equation 2 and Table 1). Red: Fit of the daily 
disc-integrated intensity excess at the Ca II K line, calculated with the relationship between Ca II K intensity excess and Bios reported 
by Harvey & White (1999). The R? between the two time series is annotated. b) The 81-day moving average. c) Zoomed inset of 


2014. 


lar intensity excess at the Ca II K line. We fit this time series to 
each of the three other disc-integrated intensity excess time se- 
ries (red, Fig. 14). The studies similar to Harvey & White (1999), 
cited in the beginning of this section, reported power-law expo- 
nents ranging from about 0.1 to 0.7. Tests showed that varying 
the power-law exponent adopted in the this calculation within 
this range makes no appreciable difference the result. 


The disc-integrated intensity excess fluctuates as network 
and faculae rotate across the solar disc due to the CLV of the 
intensity contrast/excess. At the various passbands examined in 
this study, the CLV of the intensity excess is qualitatively similar 
to the CLV of the intensity contrast (Fig. 2), which is not drawn 
here to avoid repetition. Like the intensity contrast, the intensity 
excess at the Fe I 6173 A continuum exhibits the opposite CLV 


as the Fe I 6173 A line core and 1700 A, giving the discrepant 
rotational variability in the disc-integrated intensity excess (Fig. 
14c). 

The Fe I 6173 line core and 1700 A time series are in close 
agreement with the Ca II K line time series, which is evident in 
the alignment over solar cycle (Fig. 14b) and rotation timescales 
(Fig. 14c), and in the R? (0.984 and 0.977, respectively). This 
even though the Fe I 6173 À line core and 1700 À emission do 
not conform to a similar relationship with Bios and u as Ca H 
K emission (Fig. 13). The converse is true at the Fe I 6173 A 
continuum, especially at solar rotation timescales (Fig. 14c). 

From AJ œ «By at the Ca II K line and B = By/p, it 
is straightforward to show that at a given level of B, the Ca II 
K intensity excess scales with yu. Similar to the Fe I 6173 A 
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line core and 1700 A, and opposite to the Fe I 6173 A contin- 
uum, the Ca II K intensity excess declines from disc centre to 
the limb. In terms of the time variability, the Ca II K time series 
is therefore much closer to the Fe I 6173 A line core and 1700 
A time series than the Fe I 6173 A continuum time series. We 
recall that the intensity contrast/excess CLV of the Fe I 6173 A 
continuum is opposite to that of the Fe I 6173 A line core and 
1700 A because of the different physical processes that under- 
lie the magnetic enhancement of the intensity in the lower and 
middle photosphere (Sects. 4.1.2 and 4.1.3). The results here in- 
dicate that even though the lower and middle photospheric emis- 
sion does not scale with magnetic flux density and distance from 
disc centre in the same manner as the chromospheric emission 
(Fig. 13), the time variation in the chromospheric emission from 
the network and faculae remains reasonably similar to that in 
the middle photospheric emission but not to that in the lower 
photosphere emission (Fig. 14) because of the way in which the 
CLV of the intensity contrast/excess in the various atmospheric 
regimes compare. This represents a fundamental limit to the use- 
fulness of chromospheric indices as direct indications of the ef- 
fect of network and faculae on solar irradiance. 

The disagreement between the Ca II K and Fe I 6173 À con- 
tinuum time series is caused by the opposite CLV of the inten- 
sity contrast/excess at the two passbands. In this calculation, we 
assumed the overall relationship between Ca II K intensity ex- 
cess and Bios reported by Harvey & White (1999). The authors 
also noted that when they considered the various solar surface 
magnetic feature types separately (such as enhanced network 
and active region plage), the power-law coefficients, that is, the 
constant of proportionality and the exponent, differ between the 
various feature types. The power-law exponent that Harvey & 
White (1999) reported for the various feature types they exam- 
ined ranges from 0.4 to 0.8. This means that at a given level of 
B, the Ca II K intensity excess scales with u9^ to 49? depend- 
ing on feature type, declining from disc centre to limb as in the 
overall case. We have neglected the variation in the power-law 
coefficients with feature type when we calculated the Ca II K 
time series. This likely has only very little qualitative effect on 
how the result compares to the other time series. 


4.7. Variation with time 


The 100 data days cover a nine-year period that encompasses the 
maximum of solar cycle 24,” giving us the opportunity to inves- 
tigate if the intensity contrast of network and faculae varies with 
cycle phase. For each passband and each data day, we binned the 
intensity contrast measurements by u, B, and (B), (in intervals 
of 0.25, 20 G, and 20 G, respectively) and took the binned aver- 
age. We fit a straight line to the scatter plot of the binned average 
intensity contrast and the mean over all data days. The slope of 
the straight line fit gives the scale of the intensity contrast on 
the data day relative to the mean over all data days, expressed 
in Figs. 15a. We refer to this quantity as the relative intensity 
contrast. The relative intensity contrast time series reveals that 
the intensity contrast is stable with time to about 3%. More im- 
portantly, there is no common trend between the various time 
series, nor do they indicate any clear correlation to solar cycle 
phase (by the comparison to the sunspot number and the Mg II 
index, Fig. 15b). Within the limits of our study, the intensity con- 
trast of network and faculae does not appear to vary with cycle 
phase. 


? The maximum of solar cycle 24, as indicated by the 13-month run- 
ning mean of the monthly sunspot number, was reached in April 2014. 
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Ortiz et al. (2006) repeated the analysis of Ortiz et al. (2002) 
on MDI data covering 1996 to 2001, a six-year period spanning 
the rising phase of solar cycle 23. They studied the intensity con- 
trast of network and faculae on a year-to-year basis and found it 
to be stable to about 10% over this period. Not only is this mar- 
gin of 10% non-trivial, it was not known at the time that the 
instrumental response of MDI might have changed over the ex- 
tended outages that occurred between June 1998 and February 
1999 (Ball et al. 2012). The finding here that HMI and AIA, two 
independent instruments, indicate that the intensity contrast of 
network and faculae is stable to a much smaller margin of 396 
over a longer period that encompasses almost an entire solar cy- 
cle represents a significant advance over this earlier work. 


5. Summary 


We examined the intensity contrast of network and faculae in 
the continuum and core of the Fe I 6173 À line and the AIA 
1700 À channel. We employed simultaneous full-disc observa- 
tions of the magnetic flux density and the intensity at the various 
passbands from the HMI and AIA instruments on board SDO. 
The data set comprises observations from 100 days spanning 
June 2010 to October 2018. We examined the variation in in- 
tensity contrast with distance from disc centre and magnetic flux 
density, B, the average magnetic flux density in the greater area 
about each point, called the local magnetic flux density, (Boc, 
and time. 

This is the first study of its kind to examine the variation in 
intensity contrast with local magnetic flux density over the entire 
solar disc. This was prompted by the observation that the quiet- 
Sun network appears for the same amount of magnetic flux to 
be brighter than active region faculae, a phenomena attributed to 
the fact that active regions are more magnetically crowded. We 
note that even though the heating efficiency of the magnetic flux 
is within expectation suppressed by nearby magnetic fields as a 
result of the effect of how magnetically crowded a given region 
is on flux tube sizes and how deeply embedded the magnetic 
concentrations are in intergranular lanes, the intensity contrast is 
enhanced at certain levels of B and distances from disc centre. 
Using an empirical model of the intensity contrast as a function 
of distance from disc centre, B and (B),,,, derived here from the 
intensity contrast measurements, we reproduced the difference in 
intensity contrast between the quiet-Sun network and active re- 
gion faculae. This means that the difference arises because active 
regions are more magnetically crowded alone and is not due to 
any fundamental physical difference between network and fac- 
ulae. In addition, by comparing network and faculae at similar 
distances from disc centre, B and (Boc, we found that the net- 
work and facular intensity contrast is similar to about 1%. This 
confirms that network and faculae are not fundamentally differ- 
ent from one another. The results here highlight that models of 
solar irradiance need to include the effect of nearby magnetic 
fields on the network and facular intensity. This is so far lacking 
in existing models. 

It is understood that the Fe I 6173 A line is formed largely 
in LTE, but it 1s not certain if the AIA 1700 À channel is formed 
in LTE or non-LTE. We found that we can estimate the intensity 
contrast at one passband from the measurements at the other, 
assuming that the solar atmosphere radiates like a black body. 
This means either that the AIA 1700 À channel is also formed 
in LTE, or that it is formed in non-LTE and that departures from 
LTE have little effect on the intensity contrast. If the latter is true, 
it would explain why models of solar irradiance that are based on 
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Fig. 15. a) Scale of the intensity contrast on each data day relative to the mean level over all the data days at the core (black) and 
continuum (red) of the Fe I 6173 Å line, and 1700 Å (blue). The RMS deviation from unity is indicated. b) Solar activity at the data 
days, as indicated by the sunspot number (purple) and the Mg II index (green). 


the intensity contrast of solar surface magnetic features that are 
calculated assuming LTE can reproduce the observed variabil- 
ity in spectral solar irradiance even where the LTE assumption 
breaks down. 


Chromospheric indices are employed in certain models of 
solar irradiance as direct indications of the effect of network 
and faculae on solar irradiance. The physical basis is the fact 
that chromospheric indices are largely modulated by the chro- 
mospheric network and plages that overlie the photospheric net- 
work and faculae. To investigate the validity of this approach, 
we calculated the daily disc-integrated intensity excess from net- 
work and faculae over the period of April 2010 to October 2018 
in the passbands examined here (using an empirical model based 
on the intensity contrast measurements) and in the Ca II K line 
(adopting the relationship between Ca II K intensity and mag- 
netic flux density reported by Harvey & White 1999). In terms 
of variation with time, the Ca II K time series is in close agree- 
ment with the Fe I 6173 A line core and 1700 A time series, 
but not with the Fe I 6173 A continuum time series. The rea- 
son is that the CLV of the intensity excess at the Ca II K line 
is opposite to that at the Fe I 6173 À continuum, but is qualita- 
tively similar to that at the Fe I 6173 À line core and 1700 À. 
This divergence stems from the different physical processes that 
underlie the magnetic enhancement of the intensity in the lower 
photosphere, where the Fe I 6173 À continuum is formed, and 
in the middle photosphere, where the Fe I 6173 À line core and 
the AIA 1700 À channel are formed. The results here indicate 
that chromospheric indices might be a reasonable proxy of the 
time variability in middle photospheric emission that is due to 
network and faculae, but not that in lower photosphere emission 


(Fig. 14). This represents a fundamental limit of solar irradiance 
models that employ chromospheric indices as a proxy of network 
and facular brightening. 

By comparing the intensity contrast of the network and fac- 
ulae on each data day to the average over all data days, we found 
that the intensity contrast is stable over the nine-year period ex- 
tended by the data set to about 3%. This is significantly lower 
than the margin of 10% reported by Ortiz et al. (2006). More 
importantly, there is no indication of any relationship between 
intensity contrast and solar cycle phase. This is of direct rele- 
vance to models of solar irradiance that are based on the calcu- 
lated intensity contrast of solar surface magnetic features, where 
it is assumed that the intensity contrast is invariant with time. 

The results of this study offer new and valuable insights into 
the radiant behaviour of network and faculae. These are of direct 
relevance to solar irradiance modelling. 
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